REDUCTION OF AVERAGE-TO-MINIMUM POWER RATIO 
IN COMMUNICATIONS SIGNALS 

Field of the Invention 

The present invention relates to the reduction of average-to-minimum power ratio 
in communications signals. 
State of the Art 

Many modern digital radio communication systems transmit information by 
varying both the magnitude and phase of an electromagnetic wave. The process of 
translating information into the magnitude and phase of the transmitted signal is 
typically referred to as modulation. Many different modulation techniques are used in 
communication systems. The choice of modulation technique is typically influenced by 
the computational complexity needed to generate the signal, the characteristics of the 
radio channel, and, in mobile radio applications, the need for spectral efficiency, power 
efficiency and a small form factor. Once a modulation technique has been selected for 
some specific application, it is oftentimes difficult or essentially ifripossible to change the 
modulation. For example, in a cellular radio application all users would be required to 
exchange their current mobile phone for a new phone designed to work with the new 
modulation technique. Clearly this is not practical. 

Many existing modulation formats have been designed to be transmitted with 
radios that process the signal in rectangular coordinates. The two components in the 
rectangular coordinate system are often referred to as the in-phase and quadrature (I and 
Q) components. Such a transmitter is often referred to as a quadrature modulator. (A 
distinction is drawn between a modulation, which is a mathematical description of the 



method used to translate the information into the transmitted radio signal (e.g., BPSK, 
FSK, GMSK), and a modulator, which is the physical device used to perform this 
operation.) As an alternative, the transmitter may process the signal in polar coordinates, 
in which case the signal is represented in terms of its magnitude and phase. In this case 
the transmitter is said to employ a polar modulator. A polar modulator can have several 
performance advantages over the more conventional quadrature modulator, including 
higher signal fidelity, better spectral purity, and lower dependence of device performance 
on temperature variation. 

Although a polar modulator can have practical advantages over a quadrature 
modulator, the magnitude and phase components of the signal typically have much higher 
bandwidth than the in-phase and quadrature components. This bandwidth expansion has 
implications for digital processing of the magnitude and phase, since the rate at which the 
magnitude and phase must be processed is dependent on their bandwidth. 

The rate at which the magnitude and phase varies is very much dependent on the 
modulation technique. In particular, modulation formats that lead to very small 
magnitude values (relative to the average magnitude value) generally have very large 
phase component bandwidth. In fact, if the signal magnitude goes to zero, the signal 
phase can instantaneously change by up to 180 degrees. In this case the bandwidth of the 
phase component is essentially infinite, and the signal is not amenable to transmission by 
a polar modulator. 

Many commonly employed modulation techniques do in fact lead to very small 
relative signal magnitude. To be more precise, the average-to-minimum signal magnitude 



ratio (AMR) is large. An important practical example of a modulation technique with 
large AMR is the technique employed in the UMTS 3GPP uplink (mobile-to-basestation). 

Prior work in the field may be classified into two categories: one that deals 
generally with the reduction of peak power, and another that deal specifically with "hole- 
blowing." Hole-blowing refers to the process of removing low-power events in a 
communication signal that has a time- varying envelope. This name arises in that, using 
this technique, a "hole" appears in the vector diagram of a modified signal. 

Much work has been done dealing with peak power reduction, in which the goal 
is to locally reduce signal power. By contrast, relatively little work appears to have been 
done that deals with hole-blowing (which seeks to locally increase signal power), and 
prior approaches have been found to result in less-than-desired performance. 

U.S. Patent 5,805,640 (the '640 patent) entitled "Method and apparatus for 
conditioning modulated signal for digital communications," together with U.S. Patent 
5,696,794 (the '794 patent), entitled "Method and apparatus for conditioning digitally 
modulated signals using channel symbol adjustment," both describe approaches for 
removing low magnitude (low power) events in communication signals. Both patents in 
fact refer to creating "holes" in the signal constellation. The motivation given for creating 
these holes is that certain power amplifiers, in particular LINC power amplifiers, are 
difficult to implement when the signal amplitude dynamic range is large. 

Briefly, the ' 794 patent teaches modifying the magnitude and phase of the 
symbols to be transmitted in order to maintain some minimum power. Since the symbols 
are modified before pulse shaping, the modified signal has the same spectral properties as 
the original signal. The approach used in the ' 640 patent is to add a pulse having a certain 
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magnitude and phase in between the original digital symbols before pulse shaping. 
Hence, whereas in the former patent data is processed at the symbol rate (T=l) ? in the 
latter patent, data is processed at twice the symbol rate (T=2). For brevity, these two 
methods will be referred to as the symbol rate method and the T/2 method, respectively. 
The method used to calculate the magnitude and phase of the corrective pulse(s) is nearly 
identical in both patents. 
jPJ Because the method used by both of these patents to calculate the corrective 

hi magnitude and phase is only a very rough approximation, performance is less than 

03 desired. More particularly, after processing the signal using either of these two 

O approaches, the probability of a low power event is reduced, but remains significantly 

i...J- 

;Z higher than desired. 

-11 The specific approach used in the T/2 method is to add a pulse having a 

S =3? 

Lb prescribed magnitude and phase to the signal at half-symbol timing (i.e., at t = k*T + T/2) 

before pulse shaping. The magnitude and phase of the additive pulse is designed to keep 
the signal magnitude from dropping below some desired threshold. The method does not 
allow for placement of pulses at arbitrary timing. As a result, effectiveness is decreased, 
and EVM (error vector magnitude) suffers. 

The method used in the T/2 approach to calculate the magnitude and phase and of 
the additive pulse is very restrictive in that: 

1) The signal envelope is only tested for a minimum value at half-symbol timing 
(t = i*T + T/2). 

2) The phase of the correction is not based on the signal envelope, but rather only 
on the two symbols adjoining the low-magnitude event. 
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These two restrictions can lead to errors in the magnitude and phase of the 
corrective pulses. Specifically, the true signal minimum may occur not at T/2, but at 
some slightly different time, so that error will be introduced into the magnitude of the 
corrective pulses. The validity of this assumption is very much dependent on the specific 
signal modulation and pulse-shape. For example, this may be a reasonable assumption for 
a UMTS uplink signal with one DPDCH, but is not a reasonable assumption for a UMTS 
uplink signal with two DPDCH active. The size of this magnitude error can be quite 
hi large. For example, in some cases the magnitude at T/2 is very near the desired minimum 

p;| magnitude, but the true minimum is very close to zero. In such cases the calculated 

O correction magnitude is much smaller than would be desired, which in turn results in the 

^ low-magnitude event not being removed. 

O 

Iri The signal envelope at T/2 may be greater than the desired minimum, but the 

fT signal magnitude may be below the threshold during this inter-symbol time interval, so 

that a low-magnitude event may be missed entirely. 

In any event, the correction magnitude obtained is often far from what is needed. 
The method used to calculate the corrective phase essentially assumes that the 
phase of the pulse shaped waveform at T/2 will be very close to the phase of straight line 
drawn between the adjoining symbols. This is an approximation in any case (although 
generally a reasonable one), which will introduce some error in the phase. However, this 
approximation is only valid if the origin does not lie between the previously described 
straight line and the true signal envelope. When this assumption is violated, the corrective 
phase will be shifted by approximately 180 degrees from the appropriate value. This 
typically leaves a low-magnitude event that is not corrected. 
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While the T/2 method adds pulses at half-symbol timing, the symbol rate method 
adds pulses to the two symbols that adjoin a low-magnitude event. That is, if the signal 
has a low-magnitude event at t = kT + T/2, then symbols k and (k+1) will be modified. 
Both methods calculate the phase of the additive pulses the same way, and both methods 
test for a low-magnitude event in the same way, i.e., the signal envelope at half-symbol 
timing is tested. Therefore the same sources of magnitude error and phase error 
previously noted apply equally to this method. 

The T/2 method applies the correction process repeatedly in an iterative fashion. 
That is, this method is applied iteratively "until there are no symbol interval minima less 
than the minimum threshold". 

As compared to the foregoing methods, U.S. Patent 5,727,026, "Method and 
apparatus for peak suppression using complex scaling values," addresses a distinctly 
different problem, namely reducing the Peak-to- Average power Ratio (PAR) of a 
communication signal. Large PAR is a problem for many, if not most, conventional 
power amplifiers (PA). A signal with a large PAR requires highly linear amplification, 
which in turn affects the power efficiency of the PA. Reduction is accomplished by 
adding a pulse to the original pulse-shaped waveform, with the pulse having an 
appropriate magnitude and phase such that the peak power is reduced. The pulse can be 
designed to have any desired spectral characteristics, so that the distortion can be kept in- 
band (to optimize ACPR), or allowed to leak somewhat out-of-band (to optimize EVM). 
The timing of the added pulse is dependent on the timing of the peak power, and is not 
constrained to lie at certain timing instants. 



More particularly, this peak-reduction method adds a low-bandwidth pulse to the 
original (high PAR) signal. The added pulse is 180 degrees out of phase with the signal at 
the peak magnitude, and the magnitude of the additive pulse is the difference between the 
desired peak value and the actual peak value. Because the pulse is added to the signal (a 
linear operation) the spectral properties of the additive pulse completely determine the 
effect of the peak reduction technique on the signal spectrum. The possibility is addressed 
of a peak value occurring at some time that does not correspond to a sampling instant. 

H The method emphasizes the ability to control the amount of signal splatter and/or signal 

OJ distortion. 

O U.S. Patent 6,175,551 also describes a method of PAR reduction, particularly for 

OFDM and multi-code CDMA signals where "a time-shifted and scaled reference 
IH function is subtracted from a sampled signal interval or symbol, such that each subtracted 

T 4 reference function reduces the peak power". The reference function is a windowed sine 

function in the preferred embodiment, or some other function that has "approximately the 
same bandwidth as the transmitted signal". 

Other patents of interest include U.S. Patents 5,287,387; 5,727,026; 5,930,299; 
5,621,762; 5,381,449; 6,104,761; 6,147,984; 6,128,351; 6,128,350; 6,125,103; 6,097,252; 
5,838,724; 5,835,536; 5,835,816; 5,838,724; 5,493,587; 5,384,547; 5,349,300; 5,302,914; 
5,300,894; and 4,410,955. 

What is needed, and is not believed to be found in the prior art, is a process 
whereby the AMR of a communication signal can be greatly reduced without causing 
significant degradation to the signal quality. Desirably, this process would allow the 
practical implementation of a polar modulator even for signals that have very high AMR. 
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SUMMARY OF THE INVENTION 
This invention, generally speaking, modifies pulse amplitude modulated signals to 
reduce the ratio of average power to minimum power. The signal is modified in such a 
manner that the signal quality remains acceptable. The signal quality is described in terms 
of the Power Spectral Density (PSD) and the Error Vector Magnitude (EVM). 

.y BRIEF DESCRIPTION OF THE DRAWING 

lets? 
W 

Sj The present invention may be further understood from the following description 

Si in conjunction with the appended drawing. In the drawing: 

r k ~i 

- Figure 1 is a generic block diagram for generation of a PAM signal; 

U Figure 2a is a block diagram of a PAM generator modified for AMR reduction; 

Figure 2b is a more detailed block diagram of an apparatus like that of Figure 2a; 

Figure 3 shows the impulse response of a square-root raised-cosine pulse shaping filter 

with 22% excess bandwidth; 

Figure 4 shows an I-Q plot (vector diagram) of a portion of a QPSK signal with square- 
root raised-cosine pulse shaping; 

Figure 5 shows a section of the power of the pulse-shaped QPSK signal as a function of 
time; 

Figure 6 is a histogram showing the temporal location of signal magnitude minima, 
including only those minima that are more than 12 dB below the mean power; 
Figure 7 shows phase of the QPSK signal as a function of time; 
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Figure 8 shows instantaneous frequency expressed as a multiple of the symbol rate for the 
signal of Figure 7; 

Figure 9a illustrates a situation in which a mathematical model may be used to detect the 
occurrence of a low magnitude event; 

Figure 9b is a geometric illustration of the locally linear model for the complex signal 
envelope; 

Figure 9c illustrates calculation of t_min based on the locally linear model; 
Figure 10 shows a comparison of the original signal power and the power of the signal 
after addition of a single, complex- weighted pulse designed to keep the minimum power 
above -12 dB of average power; 

Figure 1 1 is an I-Q plot of a QPSK signal that has been modified to keep the 
instantaneous power greater than -12 dB relative to the RMS power; 
Figure 12 shows the instantaneous frequency of the modified signal, expressed as a 
multiple of the symbol rate; 

Figure 1 3 shows minimum power in the modified signal as a function of error in the time 
at which the corrective pulse is added (the magnitude and phase of the corrective pulse 
are held constant); 

Figure 14 shows estimated PSD of a conventional QPSK signal and the same signal after 
application of the exact hole-blowing method; 

Figure 15a is an I-Q plot of the modified signal after demodulation, showing that some 
distortion has been introduced into the signal (the measured RMS EVM is 6.3%); 
Figure 15b illustrates the results of non-linear filtering using a root-raised cosine pulse 
that is the same as the pulse-shaping filter; 
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Figure 15c illustrates the results of non-linear filtering using a Harming window for the 
correction pulse, with the time duration equal to 1/2 the symbol duration; 
Figure 16a is a block diagram illustrating symbol-rate hole-blowing; 
Figure 16b is a more detailed block diagram illustrating symbol-rate hole-blowing; 
Figure 16c is a block diagram illustrating the iteration of symbol-rate hole-blowing; 
Figure 16d is a block diagram illustrating the iteration of symbol-rate hole-blowing; 
Pi Figure 16e is a block diagram illustrating the concatenation of one or more iterations of 

PI 

|7j symbol-rate hole-blowing followed by one or more iterations of symbol-rate hole- 

E| blowing; 

0 Figure 1 7a is a block diagram of a portion of a radio transmitter in which polar domain 

St: 

Jf; nonlinear filtering is performed; 

Ljl Figure 17b is a more detailed block diagram of an apparatus like that of Figure 17; 

1 ■■■■•£■ 
ii 'i 

jT Figure 1 8 is a waveform diagram showing the magnitude component of the polar 

coordinate signal and showing the difference of the phase component of the polar 
coordinate signal; 

Figure 19 is a diagram of the impulse response of a DZ3 pulse; 

Figure 20 is a waveform diagram showing results of nonlinear filtering of the magnitude 
component (showing the original magnitude component before polar domain nonlinear 
filtering, the magnitude component after the polar domain nonlinear filtering, the 
threshold, and the pulses added to the magnitude component); 

Figure 21 shows an example of an added pulse suitable for nonlinear filtering of the 
phase component; 
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Figure 22 is a waveform diagram showing results of nonlinear filtering of the phase- 
difference component; 

Figure 23 illustrates an alternate way for nonlinear filtering of the phase component in a 
polar coordinate system; 

Figure 24 is a block diagram of a a portion of a radio transmitter in which nonlinear 
filtering is performed first in the quadrature domain and then in the polar domain; and 
Figure 25 is a PSD showing results of the nonlinear filtering of Figure 34. 
Figure 26 is an I-Q plot for the UMTS signal constellation with one active data channel 
and a beta ratio of 7/1 5; 

Figure 27 is an I-Q plot for the UMTS signal constellation with two active data channels 
and a beta ratio of 7/1 5; 

Figure 28a illustrates the manner of finding the timings of low-magnitude events; 
Figure 28b compares probability density functions of low magnitude events using the 
exact algorithm and using the real-time approximation; 

Figure 29 is an I-Q plot illustrating a line-comparison method for vector quantization; 

Figure 30 illustrates a CORDIC-like algorithm for vector quantization. 

Figure 31 is an illustration of a known method of calculating the phase of the corrective 

pulse(s); 

Figure 32 shows Cumulative Distribution Functions (CDF) obtained with the present 
method and with known methods when the signal is pi/4 QPSK with raised cosine pulse- 
shaping and excess bandwidth of 22% (the desired minimum power is 9 dB below RMS); 



Figure 33 illustrates an example where the known symbol rate method works fairly well 
(the original signal envelope is shown, the modified envelope is shown, and the sample 
used to calculate the correction magnitude is indicated); 

Figure 34 illustrates an example where the known symbol rate method does not work 
well; 

Figure 35 illustrates an example where the known T/2 method does not work well; and 
Figure 36 shows Cumulative Distribution Functions (CDF) obtained with the present 
method and the two known methods when the signal is a 3 GPP uplink signal with one 
active DPDCH and amplitude ratio of 7/15 (the desired minimum power is 9 dB below 
RMS). 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
A polar modulator can be viewed as a combination of a phase modulator and an 
amplitude modulator. The demands placed on the phase modulator and amplitude 
modulator are directly dependent on the bandwidth of the signal's phase and magnitude 
components, respectively. The magnitude and phase bandwidth, in turn, are dependent on 
the average-to-minimum magnitude ratio (AMR) of the signal. As will be shown later, a 
signal with large AMR can have very abrupt changes in phase, which means that the 
signal phase cmponent has significant high frequency content. Furthermore, certain 
transistor technologies limit the AMR that can be achieved in a practical amplitude 
modulator. This limitation can lead to distortion of the transmitted signal if the required 
magnitude dynamic range exceeds that which can be generated by the transistor circuit. 
Thus minimization of signal AMR is highly desirable if the signal is to be transmitted 



with a polar modulator. One example of a polar modulator is described in U.S. Patent 

Application Serial No. (Dkt. ) entitled "Multi-mode 

communications transmitter," filed on even date herewith and incorporated herein by 
reference. 

The non-linear digital signal processing techniques described herein modify the 
magnitude and phase of a communication signal in order to ease the implementation of a 
polar modulator. Specifically, the magnitude of the modified signal is constrained to fall 
within a certain desired range of values. This constraint results in lower AMR compared 
to the original signal, which in turn reduces the magnitude and phase bandwidth. The cost 
of this reduction in bandwidth is lower signal quality. However, the reduction in signal 
quality is generally small, such that the final signal quality is more than adequate. 

Signal quality requirements can typically be divided into in-band and out-of-band 
requirements. Specifications that deal with in-band signal quality generally ensure that an 
intended receiver will be able to extract the message sent by the transmitter, whether that 
message be voice, video, or data. Specifications that deal with out-of-band signal quality 
generally ensure that the transmitter does not interfere excessively with receivers other 
than the intended receivers. 

The conventional in-band quality measure is the RMS error vector magnitude 
(EVM). A mathematically related measure is rho, which is the normalized cross- 
correlation coefficient between the transmitted signal and its ideal version. The EVM and 
rho relate to the ease with which an intended receiver can extract the message from the 
transmitted signal. As EVM increases above zero, or rho decreases below one, the 
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transmitted signal is increasingly distorted relative to the ideal signal This distortion 
increases the likelihood that the receiver will make errors while extracting the message. 

The conventional out-of-band quality measure is the power spectral density (PSD) 
of the transmitted signal, or some measure derived therefrom such as ACLR, ACP, etc. 
Of particular interest in relation to PSD is the degree to which the transmitted signal 
interferes with other radio channels. In a wireless communications network, interference 
with other radio channels reduces the overall capacity of the network (e.g., the number of 
simultaneous users is reduced). 

It should be clear that any means of reducing the average-to-minimum magnitude 
ratio (AMR) must create as little interference as possible (minimal degradation to out-of- 
band signal quality) while simultaneously maintaining the in-band measure of signal 
quality (i.e., EVM or rho) at an acceptable level. These considerations motivate the 
present invention, which reduces AMR while preserving out-of-band signal quality, 
which is of particular importance to operators of wireless communications networks. 

In general, AMR reduction is performed by analyzing the signal to be transmitted, 
and adding carefully formed pulses into the signal in time intervals in which the signal 
magnitude is smaller than some threshold. The details of an exemplary embodiment, 
including the signal analysis and pulse formation, are described below, starting with a 
description of a class of signals for which the invention is useful. 



Pulse- Amplitude Modulation (PAM) 

Many modern communication systems transmit digital messages using a scheme 
called pulse amplitude modulation (PAM). A PAM signal is merely a frequency- 
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upconverted sum of amplitude-scaled, phase-shifted, and time-shifted versions of a single 
pulse. The amplitude-scaling and phase-shifting of the n th time-shifted version of the 
pulse are determined by the n th component of the digital message. In the field of 
communication systems, the broad class of PAM signals includes signals commonly 
referred to as PAM, QAM, and PSK, and many variants thereof. Mathematically, a PAM 
signal x(0 at time t can be described as follows, as will be recognized by those skilled in 
the art of communications theory. The description is given in two parts, namely the 
frequency-upconversion and amplification process, and the baseband modulation process, 
as shown in Figure 1 . 

The frequency-upconversion and amplification process can be described 

mathematically as follows: 

x(t) = Re{gs(t)e JO,ct } 

where Re{} denotes the real part of its complex argument, co c ~27t/c defines the radio 
carrier frequency in radians per second and Hz, respectively, j is the imaginary square- 
root of negative unity, and g is the amplifier gain. This equation describes the frequency- 
upconversion process used to frequency-upconvert and amplify the complex baseband 
signal s(t) 9 which is also the so-called I/Q (in-phase/quadrature) representation of the 
signal. The signal s(f), which is created by the baseband modulation process, is defined 
mathematically by 

s(t) = T, a nP(t-nT) 

n 
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where p(t) is the pulse at time / and T is the symbol period (l/T is the symbol rate). For 
any time instant / at which s(t) is desired, the summation is taken over all values of n at 
which p(t-nT) is non-negligible. Also, a n is the symbol corresponding to the n th 
component of the digital message. The symbol a n can be either real or complex, and can 
be obtained from the component of the digital message by means of either a fixed 
mapping or a time-variant mapping. An example of a fixed mapping occurs for QPSK 
signals, in which the n th component of the digital message is an integer d n in the set 
{0,1,2,3}, and the mapping is given by a n =exp(jizd n /2). An example of a time-variant 
mapping occurs for 7t/4-shifted QPSK, which uses a modified QPSK mapping given by 
a„=exp(/7c(w+2rf n /4); that is, the mapping depends on the time-index n 9 not only on the 
message value d„. 

For the present invention, an important property of PAM signals is that the shape 
of the PSD (as a function off) of a PAM signal is determined exclusively by the pulse 
p{i), under the assumption that the symbol sequence a n has the same second-order 
statistical properties as white noise. This property may be appreciated by considering 
the signal s(t) as the output of a filter having impulse response p(t) and being driven by a 
sequence of impulses with weights a„. That is, the PSD S x (f) ofx(i) can be shown to be 
equal to 

= ^ - /j 2 + + /j 2 ) 

2 

where P(f) is the Fourier transform of the pulse p(t) 9 and <3 a is the mean-square value of 
the symbol sequence. 
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This important observation motivates the present method, because it suggests that 
adding extra copies of the pulse into s(t) does not alter the shape of the PSD. That is, 
nonlinear filtering performed in this manner can result in not just minor but in fact 
imperceptible changes in PSD. The adding of extra copies of the pulse into the signal can 
be used to increase the amplitude of x(f) as desired, for example when it falls below some 
threshold. Specifically, s(t) may be modified by adding additional pulses to it, to form 
new signals § (t) and x (t): 

x(t) = Re{gs(t)e J( ° ct } 

where 

S(t) = 'Za„p{t-nT)+Zb m p(t-' m l 

n m 
and the perturbation instances t m occur at points in time when it is desired to perturb the 
signal (e.g., whenever the magnitude of s(t) falls below some threshold). The 
perturbation sequence b m represents the amplitude-scaling and phase-shifting to be 
applied to the pulse centered at time t m (e.g., chosen so as to increase the magnitude of 
s(i) in the vicinity of time t m ). Like the first term in s (/), the second term in s (i) can be 
thought of as the output of a filter having impulse response p(t) and being driven by a 
sequence of impulses with weights b m . Thus, it is reasonable to expect that the PSDs of 
x (0 and x(t) will have very similar shapes (as a function of frequency/). 

With this theoretical underpinning, the present invention can be described in 
detail, in a slightly more general form than used above, as depicted in Figure 2a. The 
invention takes the signal s(t) as its input. This signal passes into an analyzer, which 
determines appropriate perturbation instances t m9 and outputs a perturbation sequence 



value b m at time instant t m . The perturbation sequence passes through a pulse-shaping 
filter with impulse response r{t\ the output of which is added to s(t) to produce s (t\ 
which in turn is passed to any appropriate means for frequency upconversion and 
amplification to the desired power. The pulse-shaping filter r(t) can be identical to the 
original pulse p(t\ as described above, or it can be different from p(t) (e.g., it may be a 
truncated version of p(t) to simplify implementation). 

A more detailed block diagram is shown in Figure 2b, illustrating a main signal 
path and a correction signal path for two signal channels, I and Q. Pulse-shaping may 
occur after pulse-shaping (sample-rate correction) or prior to pulse-shaping (symbol-rate 
correction). In the correction path, sequential values of I and Q are used to perform a 
signal minimum calculation and comparison with a desired minimum. If correction is 
required based on the comparison results, then for each channel the magnitude of the 
required correction is calculated. A pulse (which may be the same as that used for pulse- 
shaping) is scaled according to the required correction and added into the channel of the 
main path, which will have been delayed to allow time for the correction operations to be 
performed. 

The method used to determine the timing, magnitude, and phase of the correction 
pulses is dependent on the modulation format. Factors to consider include: 

1 . The duration of low-magnitude events relative to the symbol period. 

2. The timing distribution of low-magnitude events. 

If the duration of all low-magnitude events are small relative to the symbol (or chip) 
duration, then each low-magnitude event can be corrected by the addition of a single 
complex-weighted pulse, which can be identical to that used in pulse-shaping. This 



approach would be effective with, e.g., M-ary PSK modulation. The appropriate hole- 
blowing method in this case is referred to as the "exact" hole-blowing method. The exact 
hole-blowing method is described below, as well as a practical real-time hardware 
implementation of the same. Other modulation formats may lead to low-magnitude 
events of relatively long duration. This would typically be the case with, e.g., QAM and 
multi-code CDMA modulation. In such cases multiple pulses may be added, or, 
alternatively, multiple iterations of the exact hole-blowing method be used. The 
usefulness of polar-domain hole-blowing has been demonstrated to perform "final clean- 
up" of a signal produced using one of the foregoing techniques, achieving even better 
EVM performance. Because magnitude information is available explicitly in the polar- 
domain representation, hole-blowing may be performed on a sample-by-sample basis, as 
compared to (typically) a symbol-by-symbol (or chip) basis in the case of the foregoing 
techniques. 

The Exact Hole-Blowing Method 

The detailed operation of the exact hole-blowing method is now demonstrated by 
example. A QPSK signal with square-root raised-cosine pulse shaping is used for this 
example. The pulse shaping filter has 22% excess bandwidth, and is shown in Figure 3. A 
typical VQ plot of this signal is presented in Figure 4. Clearly the signal magnitude can 
become arbitrarily small. The signal power over a short period of time is shown in Figure 
5. The average power is normalized to unity (0 dB). This figure indicates that the AMR 
of this signal is at least 40 dB. In reality, the AMR for this QPSK signal is effectively 
infinite, since the signal power can become arbitrarily small. The timing of the minimum 
power is important, since the time instants at which to insert the correction pulses must be 
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determined. One would expect that the power minima would approximately occur at t = 
nT/2, where n is an integer and T is the symbol period. This is supported by Figure 5, 
which shows that the minimum power occurs very near T/2. 

To further support this hypothesis, the distribution of the power minima timing 
was examined. For this purpose, a pulse-shaped QPSK waveform with random message 
was generated, and the timing at which the power minima occurred was determined. For 
this example, it is assumed that a low power event occurs if the instantaneous signal 
power was more than 12 dB below the mean signal power. 

Figure 6 shows a histogram of the timing of the QPSK signal magnitude minima. 
These results are based on 16384 independent identically distributed symbols. Note that 
the minima are indeed closely clustered around a symbol timing of T/2, as expected. This 
is an important result, since it limits the range over which the search for a signal 
minimum must be performed. (It should be noted that the histogram shown in Figure 6 is 
only valid for this specific signal type (QPSK). Other signal types, such as high order 
QAM, may have different distributions, which must be taken into account in the search 
for local power-minima.) 

As stated previously, low-power events are associated with rapid changes in the 
signal phase. This correspondence is illustrated in Figure 7, which shows the phase of the 
QPSK signal corresponding to the power profile shown in Figure 5. It is clear that the 
phase changes rapidly near t = T/2, corresponding to the minimum power. This 
characteristic can be seen more explicitly in Figure 8, which shows the instantaneous 
frequency, defined here for the sampled data waveform as 
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where G(t) is the signal phase at time instant t and 5 is the sampling period. Figure 8 
shows that the instantaneous frequency over this interval is up to 45 times the symbol 
rate. To put this in perspective, the chip rate for the UMTS 3GPP wideband CDMA 
standard is 3.84 MHz. If the symbol rate for the QPSK signal in our example was 3.84 
MHz, the instantaneous frequency would exceed 45 X 3.84 = 172.8 MHz. Processing a 
signal with such high instantaneous frequency is not yet practical. 

It is apparent that the bandwidth of the signal phase must be reduced in order to 
enable practical implementation of a polar modulator. The most obvious approach is to 
simply low-pass filter the phase (or equivalently the phase difference). However, any 
substantial filtering of the phase-difference will lead to unacceptably large non-linear 
distortion of the signal. This distortion in turn leads to a large increase in out-of-band 
signal energy, which is typically not acceptable. Instead the recognition is made that 
rapid changes in signal phase only occur when the signal magnitude is very small. 
Therefore if the signal magnitude can be kept above some minimum value, the bandwidth 
of the signal phase will be reduced. As should be evident from the discussion of PAM 
signal spectral properties, the signal can be modified by the addition of carefully selected 
pulses without any apparent effect on the signal bandwidth. 

In order to avoid low-magnitude events, and consequently reduce the 
instantaneous frequency, a complex weighted version of the pulse-shaping filter is added 
to the signal. This complex weighted pulse is referred to as the correction pulse. The 
phase of the pulse is selected so that it combines coherently with the signal at the point of 
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minimum magnitude. The magnitude of the correction pulse is equal to the difference 
between the desired minimum magnitude and the actual minimum magnitude of the 
signal. 

Care must be taken in the calculation of the correction magnitude and phase in 
order to obtain the desired effect. Most importantly, it must be recognized that the 
minimum signal magnitude may not correspond to a sampling instant. Since the phase 
and magnitude are changing rapidly in the neighborhood of a local power minimum, 
choosing the correction phase based on a signal value that does not correspond to the 
signal minimum can lead to large phase and/or magnitude error. The likelihood of a large 
error increases as the minimum magnitude grows smaller. The likelihood of a large error 
can be lessened by using a very high sampling rate (i.e., a large number of samples per 
symbol), but this greatly, and unnecessarily, increases the computational load. 

Instead, a so-called "locally linear model" is fitted to the signal in the temporal 
neighborhood of a local power minimum, and the minimum magnitude of the model is 
solved for mathematically. The locally linear model effectively interpolates the signal so 
that the magnitude can be calculated directly at any arbitrary time instant. In this way the 
calculation of the corrective pulses are not limited to the values present in the sampled 
data waveform. This allows for highly accurate calculation of the minimum magnitude 
and, of equal importance, the required correction phase. 

Note that the complete pulse-shaped signal need not exist in order to calculate the 
correction magnitude and phase. The complex baseband signal need only be calculated at 
a small number of key sampling instants, allowing the non-linear filtering to operate on 



22 



the raw symbols before a final pulse-shaping step. This expedient can be beneficial in 
real-time implementations. 

A signal may lie outside the exclusion region defined by the desired minimum at 
two points fairly near in time, yet pass through the exclusion region in between the two 
points, as illustrated in Figure 9a. In order to find the true minimum value of the signal 
magnitude, and its corresponding phase, the following approach is preferably used. This 
approach employs a locally linear model for the complex baseband signal envelope, as 
illustrated in Figure 9b. The model approximates the complex signal envelope by a 
straight line in I/Q space in the vicinity of the local power minimum. In many cases, if 
the desired minimum power is small, and the signal modulation is not too complicated 
(e.g., PSK), this model is quite accurate. Once this model has been fitted to the signal, the 
minimum magnitude of the locally linear model can be solved for directly. 

In order to formulate the model, the complex signal must be known at no less than 
two distinct time instants. These time instants will preferably be close to the true 
minimum value of the signal, since the accuracy of the model is better in this case. These 
time instants need not correspond to sampling instants that are present in the sampled 
signal envelope. 

Denote the complex signal at time tl by s(tl). Denote the real and imaginary parts 
of s(tl) by xl and yl, respectively. Similarly, denote the real and imaginary parts of s(t2) 
by x2 and y2, respectively. The signal at these time instants are given by: 



s(/l) = 2>,p(rt-*T) 
where t2 > tl , and t2 - tl « T. Then define 

o Ax = x2 - xl 

2 Ay = y2 - 

U To the degree that the locally linear model is accurate, the complex signal at any instant 

fu in time r can be represented as 

) = + cAx + y'cAy 
Re{s(f)} = xl + cAx 
Im{s(/)}=>>l + cAy 

where c is a slope parameter. 

In order to find the minimum magnitude of this linear model, a geometric 
approach based on Figure 9b is used. The point of minimum magnitude corresponds to a 
point on the locally linear model that intersects a second line which passes through the 
origin and is orthogonal to the linear model, as illustrated in Figure 9b. This orthogonal 
line can be described parametrically (using the parameter g) by the equations: 
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x = -gAy 
y = gAx 



The point of intersection is found by setting the x-axis and y-axis components of the 
linear model and the orthogonal line equal. This yields the following set of equations 
(two equations and two unknowns): 

xl + cAx - -gAy 
y\ + cAy = gAx 

Re-arranging these equations gives: 

(gAy + xl) 



c = 



Ax 



c= (gAx-yl) 
Ay 



Setting these two expressions equal, and solving for g, gives: 



_ (ylAx - xlAy ) 
g ~ Ax 1 + Ay 2 



The minimum magnitude of the linear model is then 



( g Ax) 2 + (gAy) 2 

\g\^Ax 2 +Ay 2 

\y\Ax - xlAy\ 
4 Ax 2 + Ay 2 

The minimum magnitude as calculated from the locally linear model must be 
explicitly calculated in order to test for a low-magnitude event and to calculate the 
correction magnitude. The minimum magnitude is also used to calculate the magnitude of 
the correction pulse, which is given by: 

Pcorr ~ Pdesired — Pmm 

The correction phase is equal to the phase of the signal at minimum magnitude, as 
determined from the locally linear model. If the corrections are being performed in the 
I/Q domain, there is no need to explicitly calculate the phase 8 of the correction - only 
sin(9) and cos(9) need be calculated. Referring to Figure 9b, it may be seen that: 



2 _ 
Pmm 



Pmm 
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-gAy 



cos 6 = 



mm 



- AyjylAx - xlAy^Ax 2 + Ay : 

\y\Ax - xlAy\{Ax 2 + Ay 2 ) 

- Ay signal Ax - xlAy ) 

^Ax 2 + Ay 2 



where c/|c| = sign(c) for any scalar c. Similarly: 



. „ Ax signal Ax -xl Ay) 
sin# = , 

jAx 2 + Ay 2 



In rectangular coordinates, the in-phase and quadrature correction factors will 
respectively have the form: 



Cl= (/^desired " Anin ) COS 0 
^ = (Pde S ired-Pmin) sin6, 



Thus, it is seen that cos6 and sin6 are sufficient in order to calculate the correction 
factors, and the independent determination of signal phase 6 is not necessary. The 
modified signal in rectangular coordinates is given by: 



where t_min is the approximate time at which the local minimum occurs, and the pulse is 
assumed (without loss of generality) to be normalized such that p(0) = 1 . 

In some applications it may be desirable to obtain a precise estimate of the time 
t_min corresponding to the minimum signal power. Here, "precise" is used to mean an 
estimate that is not limited by the sample rate, and may have an arbitrary degree of 
accuracy. For example, given t_min, correction pulses may be added to the signal based 
on an interpolated or highly oversampled prototype pulse p(t). For this purpose, any of 
the three following procedures that exploit the locally linear model can be used. The first 
two approaches are based on geometric arguments. The third approach is based on direct 
minimization of the signal magnitude. 

Linear interpolation is typically employed to estimate a signal value at some 
arbitrary time t when the signal is only known at two or more discrete time indices. That 
is, linear interpolation is used to find s(t) given t. However, linear interpolation can also 
be used to find t given s(t). This property may be exploited as follows, in accordance with 
the illustration of Figure 9c. The linear model is 



s I {t) = s I (t) + c I p(t-t roiQ ) 




s(t)-s(tl) + c(Ax + jAy) 



where 



t-t\ 

c = 



tl-t\ 



For the real part of the signal, 



X 



X l + 1 i X 2 X l ) 



a t 2 -t 1 



Solving this expression for t, 



t = t l +-?—2L(t 2 -t 1 ) 



X 2 Xi 



The calculation of tjnin depends on the calculation of x_min, the real part of the 
signal at the local minimum magnitude. As described previously, the minimum 
magnitude of the signal is p min ; the real part of the signal at minimum magnitude is p min 
cosG, Therefore 
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^min _ h + 



f 



V 



Pmin COs6l - X l 
X 2 ~ X \ J 



Based on similar arguments for the imaginary part of the signal, alternatively 



1 1 j 



( 



Pmin sinO-y, 



The two expressions described above only depend on the real part or imaginary 
part of the signal In a finite precision implementation, this may be a drawback if the 
change in x or y is small. Therefore it may be better in practice to use an expression 
which is dependent on both the real and imaginary parts of the signal. Such an expression 
can be derived by direct minimization of the signal magnitude using standard 
optimization procedures. Using the linear model, the signal magnitude is 



\s{tf = (x l + cAxf + (y l + cAyf 



Taking the derivative with respect to c, 
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^-\s(tf = 2( Xl + cA x )A x + 2fo + cA>, 
oc 



Setting the derivative equal to zero and solving for c 5 



"mm 



A 2 ,+A 2 y 



Thus the final expression for t_min is 



f 



\ 



^min _ h + 
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The "exact" hole-blowing algorithm can thus be summarized as follows: 

1 . Determine the approximate timing t = tl of a potential low-magnitude event in 
the signal. 

2. In the temporal neighborhood of a potential low-magnitude event, calculate 
the pulse-shaped signal s(t) for at least two distinct time instants 

tl and t2 > tl, where t2 - 12 » T. In the case of symbol-rate hole-blowing, the 
signal s(t) is calculated based on the later-to-be-applied bandlimiting pulse 
and some number of symbols in the vicinity of tl and some number of 
symbols in the vicinity of t2. In the case of sample-rate (i.e., oversampled) 
hole-blowing, pulseshaping will have already been performed, such that s(tl) 
and s(t2) may be chosen to correspond with adjacent samples. 

3. Calculate the minimum magnitude p m i n using the "locally-linear model", 
detailed above, and t m i n? the time of p min . 

4. Compare the calculated minimum magnitude to the desired minimum 
magnitude. 

5. If the calculated minimum magnitude is less than the desired minimum 
magnitude, calculate the in-phase and quadrature correction weights Cj and 

Cq , respectively, as detailed above. 

6. Weight two copies of the pulse-shaping filter by the in-phase and quadrature 
correction values, respectively. 

7. Add these weighted copies of the pulse-shaping filter to the in-phase and 
quadrature components of the signal, referenced to t m ; n . 
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8. Translate the modified in-phase and quadrature components to magnitude and 
phase, forming the signals to be processed by a polar modulator. 

The effectiveness of the exact hole-blowing method described above is evident in 
Figure 10, compareing the instantaneous power of the original QPSK signal to the signal 
after processing with the exact hole-blowing method. The threshold for desired minimum 
power in this example was selected to be 12 dB below RMS power. It can be seen that 
the exact hole-blowing method is highly effective in keeping the signal power above the 
desired minimum value. The I-Q plot of the QPSK signal after hole-blowing is shown in 
Figure 1 L It can be seen that all traces have been pushed outside the desired limit. A 
"hole" appears in I-Q plot where none was evident before. 

The instantaneous frequency of the modified signal is shown in Figure 12. By 
comparison with Figure 8, it can be seen that the instantaneous frequency has been 
reduced from approximately 45 X the symbol rate to approximately 1 .5 X the symbol 
rate. Clearly the method has greatly reduced the instantaneous frequency. 

It should be noted that the exact hole-blowing method described here is highly 
tolerant of timing error, but relatively intolerant of magnitude and phase error (or 
equivalent errors in the correction factors ci and c Q ). Timing error refers to the time at 
which the corrective pulse is added to the original signal. This tolerance to timing error is 
due to the fact that the pulse shaping filter generally will have a fairly broad magnitude 
peak relative to the duration of the low-magnitude event. The effect of this sort of timing 
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error is shown in Figure 13, where a timing error of up to one- fourth of a symbol period 
degrades the AMR by only IdB. 

Having shown that the exact hole-blowing method is effective in removing low- 
magnitude events, the effect on in-band and out-of-band signal quality will now be 
described. The PSDs of the QPSK signal before and after hole-blowing are 
indistinguishable in the frequency domain as shown in Figure 14. The effect on EVM is 
illustrated in Figure 15a, showing an I-Q plot of the modified signal after matched- 
filtering and baud-synchronous sampling. The result shows the message that a receiver 
would obtain upon demodulation of the modified signal. It can be seen that distortion has 
been introduced to the signal, in that not all samples fall exactly on one of the four QPSK 
constellation points. The RMS EVM can be defined as 




where 62% and are, respectively, the ideal and actual PAM symbols. For the 

particular example illustrated in Figure 14, the RMS EVM was calculated using N = 
16384 symbols and found to be 6.3%. Similarly, the peak EVM is defined as follows: 
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In this example, the peak EVM is found to be about 38%. If the signal is allowed to have 
greater magnitude dynamic range, the RMS EVM will be lower. Recall, however, that 
increasing the allowed magnitude dynamic range also increases the bandwidth and peak 
value of the instantaneous frequency. Accordingly, a tradeoff can be made between the 
desired signal quality (EVM) and the instantaneous frequency requirements placed on a 
practical polar modulator. 

To assure an imperceptible effect on the signal PSD (see Figure 14), the choise of 
the correction pulse shape should substantially match the signal band limiting pulse 
shape. When the pulse oversampling rate is four or more (samples per symbol time), the 
preceeding algorithm can be simplified. Specifically, the calculation of t mm can be 
eliminated, and the correction pulse is inserted time-aligned with the existing signal 
sample s(tl). As shown in Figure 13, the maximum error in minimum magnitude is ldB 
for 4x oversampling, and decreases to 0.2dB for 8x oversampling. 

The difference in practice between an approach in which hole-blowing is 
performed taking account of t min and an approach in which the calculation of t mm can be 
eliminated is illustrated in the examples of Figures 15b and 15c, respectively. In the 
example of Figure 15b, hole-blowing was performed using a root-raised cosine pulse that 
was the same as the pulse-shaping pulse. In the example of Figure 15c, hole-blowing was 



performed using a Hanning window for the correction pulse, with a time duration equal 
to 1/2 the symbol duration. As may be seen, in Figure 15c, the original signal trajectory is 
changed as little as possible is order to skirt the region of the hole. In many instances (if 
not most), however, the calculation of t m j n can be eliminated (as in Figure 15b). 

Figure 16a shows an arrangement for symbol-rate hole-blowing. A digital 
message is applied to a main path in which pulse-shaping and upconversion are 
performed. An auxiliary path includes an analyzer block that produces a correction 
signal. An adder is provided in the main path to sum together the main signal and the 
correction signal produced by the analyzer. 

Referring to Figure 16b, a particularly advantageous embodiment for performing 
symbol-rate hole-blowing is shown. Two signal paths are provided, a main path and an 
auxiliary path. Outputs of the main signal path and the auxiliary signal path are summed 
to form the final output signal . 

The main signal path receives symbols (or chips) and performs pulse-shaping on 
those symbols (or chips), and is largely conventional. However, the main signal path 
includes a delay element used to achieve synchronization between the main signal path 
and the auxiliary signal path. 

In the auxiliary path, the symbols (or chips) are applied to a correction DSP 
(which may be realized in hardware, firmware or software). The correction DSP performs 
hole-blowing in accordance with the exact method outlined above and as a result outputs 
an auxiliary stream of symbols (or chips). These symbols (or chips) will occur at the 
same rate as the main stream of symbols (or chips) but will small in magnitude in 
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comparison, and will in fact be zero except when the signal of the main path enters or is 
near the hole. The relative timing of the main and auxiliary paths may be offset by T/2 
such that the small-valued symbols of the auxiliary path occur at half-symbol timings of 
the main signal path. 

In an exemplary embodiment, the correction DSP calculates the signal minimum 
with respect to every successive pair of symbols (or chips), by calculating what the signal 
value would be corresponding to the respective symbols (or chips) and applying the 
locally linear model. In calculating what the signal value would be corresponding to a 
symbol (or chip), the same pulse as in the main signal path and is applied to that symbol 
(or chip) and some number of previous and succeeding symbols (or chips). This use of 
the pulse to calculate what the value of the signal at a particular time will be following 
pulse-shaping is distinct from the usual pulse-shaping itself. 

After the correction symbols (or chips) of the auxiliary signal path have been 
determined, they are pulse-shaped in like manner as those of the main signal path. The 
pulse-shaped output signals of the main and auxiliary paths are then combined to form 
the final output signal. 

Since the duration of p(t) is finite (L), the signal can be evaluated ahead of time 
for the correlation between low-magnitude events and the input symbol stream (L M 
cases). The analyzer can then operate solely on the {at} to determine the c, and t,. (The 
correction symbols c ; are very small, and have negligible effect on the message stream.) 
In this variant of symbol-rate hole-blowing, the calculated symbols b m are added to the 
original message stream, and the entire new stream is bandlimited once through the filter. 



Iterative Methods 

The foregoing description has focused on symbol-rate and sample-rate hole- 
blowing techniques. Although the term "exact" has been used to refer to the these 
techniques, some inaccuracy and imprecision is inevitable. That is, the resulting signal 
may still impinge upon the desired hole. Depending on the requirements of the particular 
system, it may be necessary or desirable to further process the signal to remove these 
residual low-magnitude events. One approach it to simply specify a bigger hole than is 
actually required, allowing for some margin of error. Another approach is to perform 
hole-blowing repeatedly, or iteratively. 

In the quadrature domain, sample-rate hole-blowing may be performed one or 
multiple times (Figure 16c) and symbol-rate hole-blowing may be performed one or 
multiple times (Figure 16d). In the latter case, however, at each iteration, the symbol rate 
is doubled. For example, in a first iteration, correction symbols may be inserted at T/2. In 
a second iteration, correction symbols may be inserted at T/4 and 3T/4, etc. One or more 
iterations of symbol-rate hole-blowing may be followed by one or more iterations of 
sample-rate hole-blowing (Figure 16e). (In general, sample-rate hole-blowing cannot be 
followed by symbol-rate hole-blowing.) 

Another alternative is sample-rate hole-blowing in the polar domain. 
Polar-Domain Hole-Blowing 

The foregoing quadrature-domain hole-blowing techniques preserve ACLR very 
well at the expense of some degradation in EVM. Other techniques may exhibit a 
different tradeoff. For example, polar-domain hole-blowing preserves EVM very well at 



the expense of ACLR. Therefore, in a particular application, quadrature-domain hole- 
blowing, polar-domain hole-blowing, or a combination of both, may be applied. 

The magnitude and phase components of a polar coordinate system can be related 
to the in-phase and quadrature components of a rectangular coordinate system as 



Figure 17a is a block diagram of a portion of a radio transmitter in which polar- 
domain nonlinear filtering (i.e., "hole-blowing") is performed. The diagram shows how 
magnitude and phase as expressed in the foregoing equations are related to polar-domain 
nonlinear filtering and the polar modulator. 



In Figure 1 7a, G represents gain of the polar modulator. In operation, the digital 
message is first mapped into in-phase and quadrature components in rectangular 
coordinate system. The in-phase and quadrature components are converted into 
magnitude and phase-difference by a rectangular-to-polar converter. By knowing the 
starting point of the phase and the phase-difference in time, the corresponding phase in 
time can be calculated. This phase calculation is done at a later stage, in the polar 
modulator. Before feeding the phase-difference to a polar modulator, polar-domain 
nonlinear filtering is performed. 

A more detailed block diagram is shown in Figure 17b, illustrating a main signal 
path and a correction signal path for two signal channels, p and 0. In the correction path, 
sequential values of p are used to perform a signal minimum calculation and comparison 



p{t) = ^s 2 (t) + s 2 Q (t) 
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with a desired minimum. If correction is required based on the comparison results, then 
for each channel the magnitude of the required correction is calculated. A pulse (or pair 
of pulses in the case of the 9 channel) is scaled according to the required correction and 
added into the channel of the main path, which will have been delayed to allow time for 
the correction operations to be performed. 

An example of the magnitude and phase components (phase-difference) in polar 
coordinates is shown in Figure 18. As can be seen from the plot, when the magnitude 
dips, there is a corresponding spike (positive or negative) in the phase-difference 
component. 

The spike in the phase-difference component suggests that there is a rapid phase 
change in the signal, since the phase and phase-difference have the following 
relationship: 



The dip in magnitude and the rapid phase change are highly undesirable as both 
peaks expand the bandwidth of each polar signal component. The goal of the polar 
domain nonlinear filtering is to reduce the dynamic range of the amplitude swing as well 
as to reduce the instantaneous phase change. The dynamic range of the amplitude 
component and the maximum phase change a polar modulator can handle is limited by 
the hardware. Polar domain nonlinear filtering modifies the signal before the signal 
being processed by the polar modulator. This pre-processing is to ensure that the signal 
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dynamic range will not exceed the limit of the implemented hardware capability of the 
polar modulator, and therefore avoids unwanted signal distortion. 

Nonlinear filtering in a polar coordinate system is somewhat more complex than 
nonlinear filtering in a rectangular coordinate system. Both components (magnitude and 
phase) in the polar coordinate system have to be handled with care to prevent severe 
signal degradation. 

The following description sets forth how polar-domain nonlinear filtering is 
performed in an exemplary embodiment. Polar-domain nonlinear filtering is composed of 
two parts. The first part is nonlinear filtering of magnitude component, and the second 
part is nonlinear filtering of the phase component (phase-difference). These two parts will 
each be described in turn. 

Nonlinear filtering of magnitude component 

If the magnitude dynamic range exceeds the capability of the polar modulator, the 
output signal will be clipped. This clipping will result in spectral re-growth and therefore 
greatly increase the adjacent channel leakage ratio (ACLR). One method that can be 
applied to reduce the dynamic range of the magnitude component is hole-blowing (or 
nonlinear filtering). 

The purpose of this nonlinear filtering of the magnitude component is to remove 
the low magnitude events from the input magnitude p(t), and therefore reduce the 
dynamic range of the magnitude swing. Assume a threshold (THmag) for the minimum 
magnitude; by observing the signal p(t), it is possible to obtain the time intervals where 
the signal falls below the threshold. Assume there are N time intervals where the signal 



falls below the threshold and that the magnitude minimums for each of these time 
intervals happen at ti y ... , tN, respectively. Therefore, the nonlinear filtering of 
magnitude component can be expressed as: 

N 
n=l 

where bn and p n (t) represent the inserted magnitude and pulse for the n th time interval. 
The nonlinear filtered signal is composed of the original signal and the inserted pulses. 
The magnitude of the inserted pulse is given as: 

b n =TH mag -p{t n ) 

The inserted pulse p n (t) has to be carefully chosen so that the signal degradation 
with respect to ACLR is minimized. It is desirable for the pulse function to have smooth 
leading and trailing transitions. A suitable pulse (DZ3) is based on McCune's paper 
entitled "The Derivative-Zeroed Pulse Function Family CIPIC Report #97-3, 
University of California, Davis, California, June 29, 1997. The impulse response of the 
DZ3 pulse is plotted in Figure 19. 

Results of an example for nonlinear filtering of the magnitude component are 
shown in Figure 20. As can be seen from the plot, signal intervals where the magnitude 
dropped below the threshold are compensated by the inserted pulses. As a result, the 
dynamic range of the magnitude swing is reduced. 
Nonlinear filtering of phase component 
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In an exemplary embodiment, the input to the polar modulator is phase- 
difference, not phase. A voltage-controlled oscillator (VCO) in the polar modulator 
integrates the phase-difference and produces the phase component at the output of the 
polar modulator. The phase-difference is directly related to how fast the VCO has to 
integrate. If the phase-difference exceeds the capability of the VCO, the output signal 
phase will lag (or lead) the actual signal phase. As a result, phase jitter will occur if the 
VCO is constantly unable to keep up with the actual signal phase. This phase jitter will 
result in constellation rotation, and may therefore severely degrade EVM. 

The purpose of the nonlinear filtering of phase component is to suppress large 
(positive or negative) phase-difference events so that phase accumulation error will not 
occur. It is important to know that the output of the VCO is an accumulation of the input 
(phase-difference). Therefore, any additional processing to the phase-difference has to 
ensure that the phase error will not accumulate. The nonlinear filtering of the phase 
component described presently carefully modifies the phase-difference in a way that the 
accumulated phase might deviate from the original phase trajectory from time to time. 
However, after a certain time interval, it will always merge back to the original phase 
path. 

Nonlinear filtering of the phase component is done by first finding the location 
where the absolute value of the phase-difference is beyond the integral capability of the 
VCO. Assume there are a total of M events where the absolute values of the phase- 
difference are beyond the capability of the VCO, and that the peak absolute value for 
each event happens at tm. Nonlinear filtering of the phase-difference component may 
then be described by the following expression: 
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M 

e(t)=e(t)+^c mPp ^,(t-t H ) 

m=l 

where p Pl m(t) is the pulse inserted to the phase-difference component at time tm 9 and Cm is 
the corresponding magnitude given by 

c m =m p -e(t m ) 

where TH P is the threshold for the phase-difference. 
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The inserted pulse should satisfy the following equation 



The result of inserting a pulse is basically the same as altering the phase 
trajectory. However, by inserting pulses that satisfy the above equation, the modified 
phase trajectory will eventually merge back to the original phase trajectory. This can be 
seen by the following equation: 



m=l 

The second term on the right-hand side will eventually disappear. Therefore, the 
modified phase trajectory will merge back to the original phase trajectory. 

An example of an added pulse p P (t) suitable for nonlinear filtering of the phase 
component is plotted in Figure 21. Again, it is very important for the pulse function to 
have smooth leading and trailing transitions. The pulse used for nonlinear filtering at the 
phase-difference path is composed of two pulses. These two pulses have the same area 
but different polarity and durations. Therefore, the integration of the combined pulse 
with respect to time is zero. 

Results of an example of nonlinear filtering of the phase-difference component 
are shown in Figure 22. As can be seen from the plot, the modified phase trajectory has a 



smoother trajectory than the original phase trajectory. The modified phase trajectory 
requires less bandwidth from the VCO than the original phase trajectory. 

Polar-domain nonlinear filtering is composed of the filtering of magnitude and 
phase components. If the nonlinear filtering is done jointly (both magnitude and phase), 
better spectral roll-off can be achieved. However, each nonlinear filtering operation can 
also be done independently. 

Instead of nonlinear filtering of the phase-difference component, direct nonlinear 
filtering of the phase component can also be implemented. Figure 23 shows an example 
in which the original signal has a steep increase in phase from time ti to ti. One way of 
reducing the phase change is by interpolating between the straight line v(t) and the 
original phase path. The interpolation can be expressed as: 

§(t) = w(t)0(t) + (1 - w(0)v(0 

where w(t) is a weighting factor. The weighting factor can be a constant or may vary in 
accordance with DZ3, a Gaussian function, or the like. Again, it is important for the 
weighting function to have smooth leading and trailing edges. 

Polar-domain nonlinear filtering can be used in concatenation with quadrature- 
domain nonlinear filtering. Such an approach may involve less computational 
complexity than iterative quadrature-domain nonlinear filtering. 

If only one iteration of quadrature-domain nonlinear filtering is performed, low- 
magnitude, high-phase-variation events may still occur, but with low probability. These 
low probability events can degrade signal quality if not properly managed. To eliminate 



these low-probability events, the foregoing polar-domain nonlinear filtering method may 
be used following a single iteration of quadrature-domain nonlinear filtering, maintaining 
low EVM and low ACLR. Figure 24 shows the block diagram for this concatenated 
system. Figure 25 shows the PSD of the output signal processed by the quadrature-polar 
nonlinear filtering module. Curve A represents the signal with eight iterations of 
quadrature-domain non-linear filtering. Curve B represents the signal with one iteration 
of quadrature-domain nonlinear filtering followed by one iteration of polar-domain 
nonlinear filtering. Spectral re-growth is below -60dB. In addition, good spectral roll-off 
is achieved. 

Reduction of Computational Load for Real-Time Applications 

In a practical polar modulator, the hole-blowing algorithm must be implemented 
in real-time using digital hardware and/or software. Real-time implementation of the 
exact hole-blowing method presents several challenges. The particular challenges faced 
are, to some extent, dependent on the overall architecture selected. There are (at least) 
two alternative architectures for implementing the exact hole-blowing algorithm. In the 
first architecture, referred to as symbol-rate hole-blowing, corrective impulses are 
calculated and added to the data stream with appropriate timing, after which pulse 
shaping is performed. In the second architecture, referred to as sample-rate (or 
oversampled) hole-blowing, the full pulse-shaped signal is calculated, after which 
weighted pulses are added to the pulse-shaped signal. 

Because of its computational complexity, the exact hole-blowing algorithm 
described earlier is difficult to implement using digital hardware. In the exact hole- 
blowing algorithm, arithmetic division, square, and square-root operations are required. 



These arithmetic operations greatly increase the complexity of the digital hardware and 
should be avoided if possible. Reduction in computational complexity of the 
implemented algorithm leads directly to a reduction in hardware complexity. The real- 
time hole-blowing algorithm described herein requires no division, square, or square-root 
operations. 

Although a real-time hole-blowing algorithm can be implemented in either 
symbol-rate or sample-rate form, it is generally more desirable to implement symbol-rate 
hole-blowing. In particular, symbol-rate hole-blowing algorithm allows the digital 
hardware to be operated at a slower clock speed and does not require as much memory to 
perform pulse insertion. A symbol-rate real-time hole-blowing algorithm will therefore 
be described. 

Generally speaking, real time hole-blowing algorithm is very similar to the exact 
hole-blowing algorithm. However, some assumptions and approximations are made in 
real time hole-blowing algorithm to simplify implementation. In the real-time hole- 
blowing algorithm, assumptions are made regarding where magnitude minimums are 
most likely to occur based on the structure of the signal constellation(s). The assumption 
allows for determination of where to insert the pulse, if needed, without calculating the 
entire signal magnitude. In addition, the arithmetic operations in exact hole-blowing 

I 2 2 

algorithm are greatly simplified by normalizing the expression <\j Ax + Ay to one. 
This normalization eliminates the need for division, square, and square-root operations. 

The differences between the real-time and exact hole-blowing algorithms may be 
appreciated by considering answers to the following questions: 

1 . What are the potential timings for the low-magnitude events? 
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2. What are the values of magnitude minimums, and how can they be calculated 
efficiently and accurately? 

3. What are the in-phase and quadrature correction weights of the magnitude 
minimum if a signal falls below the prescribed threshold? 

These questions will be addressed using the UMTS signal as an example. 
L What are the potential timings for the low-magnitude events? 

If the timings of the magnitude minimums can be estimated based on the input 
data bits, then it is not necessary to calculate the whole waveform in order to obtain the 
magnitude minimums. This shortcut affords a great saving in computation. The locally 
linear model can be used to calculate the minimums with great accuracy if the 
approximate timing for low-magnitude events is known. 

As previously described, the minimum magnitudes for a bandlimited QPSK signal 
usually occur close to half-symbol instance (nT+T/2). This assumption is supported by 
the histogram in Figure 6. This assumption can also be applied to higher-order pulse- 
shaped PSK signals whose constellation points lie on the same circle. A good example 
for this type of signal is the UMTS signal with one active data channel. Figure 26 shows 
an example of the UMTS signal constellations with one active data channel. As can be 
seen from the plot, all constellation points lie on the same circle. If a histogram is 
constructed for this particular signal, it may be seen that the timings for the low- 
magnitude events are very likely to happen at every half-symbol instance. Therefore, for 
this particular signal, it is assumed that the timings for the low-magnitude events will 
occur at time nT+T/2, where n is an integer. 



The above assumption will not be true in the case of a more complex signal 
constellation. Figure 27 shows an example of the UMTS signal constellations with two 
active data channels. As can be seen from the plot, not all the constellation points lie on 
the same circle. The algorithm used to find the magnitude minimums for a more complex 
signal constellation is illustrated in Figure 28a, and described as follows: 

1 . If the signal transitions from constellation point PI at time nT and ends at 
point P2 at time (n+l)T, a straight line is drawn connecting the two points. 
This line is denoted as the first line. 

2. A second line is drawn perpendicular to the line that connects PI and P2. The 
second line includes the origin and intersects the first line at point M. 

3. The second line divides the first line into two sections whose length are 
proportional to Dl and D2, respectively. The timing for the low-magnitude 
event is approximately nT+Dl/(Dl+D2)T. 

4. If the second line and the first line do not intersect, there is no need for pulse 
insertion. 

5. Having determined where the magnitude minimums are most likely to occur, 
the locally linear model is then used to calculate the local magnitude 
minimums. 

With the above algorithm, a look-up table can be built for different signal 
constellations. The size of the look-up table can be reduced if the constellation points are 
symmetrical to both x and y axis. A special case of the foregoing algorithm occurs where 
Dl always equal to D2, corresponding to the UMTS signal with one active data channel. 
Therefore, the magnitude minimums for one active data channel always happens close to 
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time nT+T/2. The above algorithm for finding the approximate locations of magnitude 
minimums can also be generalized to signals with more complex constellations. 

Figure 28b shows the probability density function, by sample interval (assuming 



both the exact method described earlier and the real-time approximation just described. 
Note the close correspondence between the two functions. 

2. What are the values of the magnitude minimums. and how can they be calculated 
efficiently and accurately? 

With the knowledge of the timing where the magnitude minimums are likely to 
occur, the locally linear model is then used to calculate the magnitude minimums. A 
determination must be made whether the magnitude of the signal falls below the 
prescribed threshold. If the magnitude falls below the prescribed threshold, in-phase and 

quadrature correction weights, Cj and Cq , for the inserted pulse have to be calculated. 

From the exact hole-blowing algorithm, the magnitude minimum can be 
calculated by the following equation: 



15 samples per symbol), of minimum magnitude events where correction is required, for 




where xl and yl are the in-phase and quadrature samples of the signal at time tl, and 
tl is the time where the magnitude minimums are most likely to occur. 



The equation above involves one division, one square-root operation, two square 
operations, as well as several other operations ~ multiplication, addition and subtraction. 
The computational complexity can be reduced by normalizing the denominator of the 
above equation. 

Let the vector (Ax, Ay) be expressed as: 



Ax = /^cos(#J 
Ay = p x} sm(6J 



where p^ and 0^ are the magnitude and phase of the vector (Ax, Ay 7 ). Then, the 
minimum magnitude of the signal can be expressed as 

Pm'm ~ 



^Ax 2 + Ay : 



XV 



jlcos(^)-xlsin(^)| 



Pxy V COs2 ( 0 xy ) + sin2 Wxy > 



With this normalization, no division, square, or square-root operation is required. 
In addition, it is not necessary to know p^ . However, the values of sin(0^ ; ) and 

COSiPxy ) are needed. Therefore, given a vector (Ax, Ay), a way is needed to obtain 

(sin(#^ ), COS^?^ )) efficiently. 
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Two ways will be described to evaluate sin^^) and COS(^) with 

relatively low hardware complexity. The first one uses a line-comparison method, and 
the second one uses the CORDIC (Coordinate Rotation for Digital Computer) algorithm . 

First, consider the approximation of sin(#^) and COS^^) using a line- 
comparison method. As can be seen from Figure 29, lines with a function of y = Mx 
partition the first quadrant into several sub-sections. By comparing the point 
(| Ax |, | Ay |) against the lines y = Mx , the sub-section that the point 
(| Ax |, | Ay |) falls into can be determined. Any point within a certain sub-section is 
represented by a pre-normalized point (sin(^ ), COS(#^ ) ), where i denotes the 

section the point (| Ax |, | Ay |) belongs to. 

The details of the algorithm are described as fellows: 

1 . First, convert the vector (Ax, Ay) into the first quadrant (| Ax \, \ Ay |). 

2. Compare | Ay | with M \ Ax | where M is positive. 

3. If | Ay | is greater than M | Ax |, then (| Ax |, | Ay |) is located at the left 
side of the line y ~ Mx . 

4. If | Ay | is smaller than M \ Ax \, then (| Ax |, | Ay |) is located at the 

right side of the line y = Mx . 

5. Based on the above comparisons with different lines, the sub-section the point 
(| Ax |, | Ay |) belongs to can be located. Assuming the point belongs to 
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section i , then the vector (sin^^ ), COsC^ )) can be approximated by 

the vector (sign(Ax) * sin(<9^ ), sign(Ay) * cos(<9^ )). 

A table is used to store the pre-calculated values of (sin(^),cos(^)). if 

a total of W lines are used for comparison, there will be (W+l) sub-sections in the first 
quadrant. As a result, the whole vector plane is being divided into 4*(W+1) sub-sections. 
The second method used to approximate (sin(^ jqy ), COS(^ )) is a CORDIC- 

like algorithm. This method is similar to the line-comparison method. However, the 
CORDIC-like algorithm partitions the vector plane more equally. A more detailed 
description of the algorithm is presented as follows and is illustrated in Figure 30: 

1 . First, convert the vector (Ax, Ay) into the first quadrant (| Ax |, | Ay |). 

2. Second, rotate the vector (| Ax |, | Ay |) clock-wise with an angle 
0 O = tan" 1 (1) . The vector after the angle rotation is denoted as 

(| Ar| 0 ,|Ay | 0 ). 

3. Let/ = 1. 

4. If | Ay is greater than 0, rotate the vector (| Ax , \ Ay ) 
clockwise with an angle 0 Z = tan -1 (2~ l ) . Otherwise, rotate the vector 
counter-clockwise with an angle of 0; . The vector after the angle rotation is 

(|A*WAj;| z .). 

5. Let i = / + l. 



6. Repeat steps 4 and 5 if needed. 

7. Assume a total of K vector rotations are performed. This algorithm 
partitions the first quadrant into 2 A k sub-sections. The sign of 

j Ay |o , | Ay |j 9 • • • , j Ay can be used to determine which of the 2 A k 

sub-sections the vector (| Ax | ? | Ay |) belongs to. 

8. A table filled with pre-quantized and normalized values is then used to 
approximate the vector (| Ax \ 9 \ Ay |). 

9. If the look-up table gives a vector of (sin(#^, ) ? COSf^ ) ) for vector 

(| Ax |,| Ay |), then (sin(0^ ), COS(6 xy )) can be approximated by the 
{sign(Ax) * sin(^ ), sign(Ay) * cos(^ )). 



vector 



The vector rotation in the CORDIC algorithm is carefully done so that the vector 
rotation is achieved by arithmetic shifts only. This leads to a very efficient structure. 
The accuracy of the approximation can be improved by going through more CORDIC 
iterations. If a total of two vector rotations are performed, the resulting partition of the 
first quadrant is similar as in the line-comparison method shown in Figure 29. 



I 2 2 

The normalization of the expression ^ Ax + Ay and the efficient algorithms 
to approximate (sill^^ ), COS^^ )) greatly reduce the computational complexity of 
the locally minimum method. Using these methods, p m ^ Q may be readily evaluated. If 
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the magnitude minimum P m [ n falls below the threshold P desired > a P u l se insertion 
will be performed based on the hole-blowing algorithm. This leads to the third question: 
3, What are the in-phase and Quadrature correction weights if the magnitude minimum of 
the signal falls below the prescribed threshold? 

Similar techniques can be applied to the calculation of sill(#) and COs(6?) by 
replacing Ax with Ax — p^ 00^(9^) and Ay with Ay = p^ sin(# xj; ) , as 
follows: 



sin# = 



Axsign(jlAx - xlAy) 



^Ax 2 +Ay 2 
Pxv cos(0 w ) signal Ax - xlAy) 




= cos(^ ) signal Ax - xlAy) 
cos 0 = - sin(6> ) signal Ax - xl Ay ) 



With the above simplification, the in-phase and quadrature correction weights. 



Cj and Cq , can be calculated easily with simple subtraction and multiplication: 
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Summary of the real time hole-blowing algorithm 

To facilitate real-time implementation, a novel approach is taken to estimate the 
timings for the low-magnitude events. Also, the locally-minimum method is greatly 

2 2~ 
Ax +Ay . As a result, no division, 

square operation, nor square-root operation are required. 
M Two methods are being proposed to evaluate (sin(0^),COS(0^)) given the 

; j: 

bj vector (Ax, Ay). The first one is the line-comparison method, and the second one is the 

S3 CORDIC-like algorithm. The implementation complexities of the above two methods are 

u generally low since only arithmetic shifts and comparisons are needed. 

ii ;: 
p~ 
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APPENDIX 1 

A1.0 Analysis of the O'Dea Hole-Blowing Methods 

In this Appendix, a detailed analysis is presented of the hole-blowing methods proposed 
in U.S. Patents 5,696,794 and 5,805,640. These two patents are very similar and are 
therefore treated simultaneously. The main difference is that the former patent 
(5,696,794) modifies the symbols to be transmitted, while the latter patent (5,805,640) 



adds pulses at T/2 symbol timing instants. For brevity, patent 5,696,794 will be referred 
to as the symbol rate method, and patent 5,805,640 will be referred to as the T/2 method. 
An overview of both methods is first presented, followed by an examination of 
performance with two different signal modulations. The first test signal is 7i/4 QPSK with 
zero-ISI raised cosine pulse shaping. This is the modulation employed in the two patents. 
The second test signal is a UMTS 3 GPP uplink signal with one active DPDCH and a 
DPDCH/DPCCH amplitude ratio of 7/15. UMTS uses square-root raised-cosine pulse 
shaping with 0.22 rolloff characteristic. 

Al.l Overview of the Q'Dea Hole-Blowing Algorithms 

The term "half-symbol timing" is defined to be those time instants that are exactly 
halfway between symbol times. That is, if the PAM signal is modeled as 




58 



where T is the symbol period and p(t) is the pulse shape, then the half-symbol times 
correspond to t = kT + T/2 where k is an integer. For clarity of presentation, it will be 
assumed that the maximum value of p(t) has been normalized to unity. 

Both methods test for the existence of undesirable low-power events by measuring the 
signal magnitude at the half-symbol time instants, and comparing this value to some 



Both patents use the same method to calculate the phase of the corrective pulse(s). 
Assume that the low-magnitude event occurs between symbols k and &+1 . First, 
determine the so-called phase rotation 0 rot , which is simply the change in phase in the 
transition from symbol & to symbol £+1, as illustrated in Figure 31. The corrective phase 
6adj is given by 



where Ok is the phase of the kth symbol. A vector with phase equal to the adjustment 
phase is orthogonal to a straight line drawn from symbol k to symbol £+1, as illustrated in 
Figure 31. Note that since there are a finite number of possible phase rotations, there are 
also a finite number of possible phase adjustments, so that explicit calculation of the 
rotation phase is not necessary. 



desired minimum magnitude mag_d: 



? 




9adj ~®k 
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If the symbol rate approach (5,696,794) is used, the two symbols adjoining the low- 
magnitude event (i.e., symbols k and £+1) are modified by the addition of a complex 
scalar. The magnitude of this complex scalar is given by 



where p mid is the amplitude of the pulse shaping filter at t = T/2. (The rationale for 
calculating the magnitude of the correction in this manner, however, is unclear.) The 
complex scalar adjustment is then 



Note that both symbols are modified in the same manner. 

Thus, a noise component is added to the signal by intentional relocation of the 
information symbols {a k }. This can "confuse" any equalizer at the receiver, which will 
expect that any signal distortion is due to the channel. 

If the T/2 approach (5,805,640) is used, a complex scalar is added to the symbol stream, 
before pulse-shaping, at the appropriate half-symbol time instants. When a low 
magnitude event is detected at t = kT + T/2, a complex symbol with magnitude 





and the resulting modified symbols are given by 



a k =a k~ ]r C adj 
a k+l = a k+\ + C adj 



m = (magd - mag_s) 

and phase equal to e a dj (given above) is added at t = kT+T/2. 

This restriction to T/2 insertion timing limits this method to circular signal constellations. 

One additional difference between the symbol rate method and the T/2 method is that the 
symbol rate method is intended to be applied iteratively until the signal magnitude does 
not drop below some threshold. There is no mention of an iterative process in the T/2 
patent. 

Al .2 Performance with id A QPSK 

Consider now the performance of the known hole-blowing algorithms relative to the 
disclosed "exact" hole-blowing method when the targeted signal is k/4 QPSK. It is 
interesting to note that this signal has a "hole" in its constellation if the rolloff is high, 
e.g., a=0.5. A rolloff of 0.22 was chosen so that the signal would not have a pre-existing 
hole. 

Figure 32 shows the CDF obtained from the disclosed method and the two known 
methods with the aforementioned p/4 QPSK signal. The desired minimum power level 
was set at 9 dB below RMS. These simulation results are based on 16384 symbols with 
32 samples/symbol. The figure clearly shows that the exact method is much more 
effective than either of the known methods. It is also evident that both known methods 
perform similarly. This is not surprising given the similarity of the two approaches. 



Some explanation of the performance of the known methods as shown in Figure 32 is in 
order. It was noted earlier that there are sources of error possible in calculation of both 
the magnitude and phase of the correction pulses. Figure 33 shows an example where the 
prior-art symbol rate method works fairly well. The signal envelope is not pushed 
completely out of the desired hole, but the method is performing more-or-less as 
intended. In contrast, Figure 34 shows an example where the prior-art symbol rate 
method does not perform well This example illustrates error in both the magnitude and 
phase of the corrective pulses. In this example, the trace passes on the "wrong side" of 
the origin relative to the assumptions made in calculating the correction phase. Thus the 
trace is pushed in the wrong direction. Furthermore, it is evident from this example that 
the magnitude at T/2 is not the minimum magnitude, so that even if the phase had been 
calculated correctly, the signal would not have been pushed far enough. 

Figure 35 shows an example where the prior-art T/2 method does not perform well. (The 
segment shown in Figure 35 is the same segment of the signal shown in Figure 34, where 
the symbol rate method did not perform well.) Comparing Figure 34 and Figure 35, it is 
clear that the two methods yield nearly identical traces. It can be seen that the symbol rate 
method only alters the symbols adjacent to the low-magnitude event, while the T/2 
method affects more symbols. 
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A1.3 Performance with a 3 GPP Uplink Signal 

Consider now the performance of the known methods with a more realistic signal, that 
being a 3 GPP uplink signal with one active DPDCH and an amplitude ratio of 7/1 5. 
Figure 36 shows the CDF's obtained with the disclosed exact method and with the known 
hole-blowing methods when applied to one frame (38400 chips) of the signal with 32 
samples/chip. It is clear that the exact correction method greatly outperforms the known 
methods, and that the known methods perform comparably. 



